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Abstract 

Over the last three years, several satellite and balloon observatories have suggested intriguing features in the cosmic 
ray lepton spectra. Most notably, the PAMELA satellite has suggested an "anomalous" rise with energy of the cosmic 
ray positron fraction. In this article, we summarize the global picture emerging from the data and recapitulate the 
main features of different types of explanations proposed. The perspectives in testing different scenarios as well as 
inferring some astrophysical diagnostics from current/near future experiments are also discussed. 
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1. Introduction 

Probes of very different nature, both cosmological 
and astrophysical, suggest that the Universe is CP- 
asymmetric, containing negligible amounts of antimat- 
ter. Concerning our Galactic environment, this is con- 
firmed by the strong dominance of the matter with re- 
spect to the antimatter component in cosmic rays (CRs). 
The small amounts of positrons and antiprotons de- 
tected in CRs are attributed to byproducts of collisions 
of CRs in the rarefied interstellar medium (ISM), a pro- 
cess whose probability can be inferred from the gram- 
mage deduced by other secondary/primary ratios like 
the boron-to-carbon ratio, B/C. In turn, this permits to 
use these species as additional handles on CR propaga- 
tion and interactions in the ISM as well as to perform 
consistency checks of the models. 

Turning the argument around, under the assumption 
of purely secondary production, the positron flux can be 
computed in a relatively robust way, of course subject to 
some uncertainties coming from primary fluxes, cross 
sections and propagation parameters, see e.g. JH. Devi- 
ations from these expectations have even been proposed 
as indirect signatures of dark matter (DM) in the Galac- 
tic halo (see the pioneering article (2[] and refs. to it): In 
most popular models where DM is made of (meta)stable 
weakly interacting massive particles (WIMPs), their an- 
nihilations or decays in the Galactic halo produce GeV- 
TeV energy cosmic rays equally abundant in particles 
and antiparticles. Since antiprotons constitute only 10~ 4 



of the proton flux and positrons only ~ (3(10%) of the 
electron flux, the lower background makes antimatter 
the preferred channel for searches of WIMP signatures 
in the flux of charged cosmic rays. 

It is thus not surprising that a great excitation has 
followed the release of new data on the positron frac- 
tion e + /(e + + O by the PAMELA collaboration 0,0], 
as well as the important complementary information on 

e + + e~ flBHll&H] and e ~ s P ectra id- In Par- 
ticular, the positron fraction in the cosmic ray spectrum 
measured by PAMELA appears to begin climbing quite 
rapidly between ~7 GeV till at least -100 GeV. Al- 
though previous experiments such as HEAT ill 311 and 
AMS-01 II 1 411 were already hinting to a possible devia- 
tion from basic expectations, the PAMELA data make 
this evidence striking thanks to a much higher statistics 
and extend it over a wider energy range. A recent, in- 
dependent confirmation of this behaviour has been pub- 
lished by Fermi lfl5ll . 

The concern was raised (see e.g. fljlo that this signal 
could be mimicked by misidentified protons in the in- 
strument, if a rejection power no better than ~ <9(10~ 4 ) 
was achieved by PAMELA. Since then, the collabora- 
tion has extensively argued (see [4]) that their rejection 
power — both based on test beam data and on in-flight 
data — can be estimated at the 10~ 5 level or lower, thus 
strongly disfavouring proton contamination as a plausi- 
ble cause. In the following, we shall assume that the 
data reflect a real feature. This will be anyway indepen- 
dently confirmed by AMS-02 11711 . which was launched 



Preprint submitted to Astroparticle Physics 



June 7, 2012 



on May 16, 201 1 and is successfully operating onboard 
of the International Space Station. Note that a more sub- 
tle artifact signal may follow if the detected particles 
are truly positrons, but rather deriving from locally pro- 
duced secondary particles, for example in the pressur- 
ized container of the PAMELA detector. Judging from 



the simulations reported in 111 811 . it appears that the num- 
ber of such events is negligible. 

Here we have no ambition of reviewing in depth 
the observational situation of leptonic CR data (a re- 
cent review covering these aspects has been presented 
in fl^l .) In particular, for a detailed discussion of the 
error budget we refer to the original publications. Yet, 
since we shall refer to recent observations extensively, 
in Fig. Q] we report the PAMELA and Fermi positron 
fraction data together with a couple of previous deter- 
minations at low energy, while in Fig. [2] we report the 
electron+positron flux as published by several experi- 
ments since 2008 together with the electron-only flux 
recently presented by PAMELA and the electron-only 
and positron only fluxes recently presented by Fermi. 
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Figure 1: The positron fraction vs. energy measured by PAMELA 
2008 @1 (statistical errors only), PAMELA 2010 01 (systematical 
and statistical errors summed in quadrature), and Fermi 2011 [15], 
compared with previous results from HEAT (combined data) [Till and 
AMS-01 Q]||. 

This review article is structured as follows: in Sec. [2] 
we sketch why these data appear to require a source 
of primary positrons (or, most likely, of electrons and 
positrons). In Sec. [3] we briefly summarize why DM is 
unlikely to be involved in the explanation of the phe- 
nomenon. In Sec. |4] we outline the basic arguments 
in favour of pulsars as sources, as well as some open 
problems in the astrophysics of these objects. In Sec. [5] 
we describe the main alternative, namely (hadronic) 
production and acceleration in Supernova Remnants 
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Figure 2: Electron plus positron fluxes published since 2008 by PPB- 
BETS | 5] (statistical errors only), ATIC [6] (statistical errors only), 
HESS (71 9] (statistical errors only) Fermi |8, 10], as well as the e~ 
spectra recently published by PAMELA 1 1 1 ] and the e~ and e + spec- 
tra recently presented by Fermi [15]. For Fermi and PAMELA data, 
systematics and statistical errors are summed in quadrature. 



(SNRs), which are considered the standard sources of 
hadronic CRs 11211 . In Sec. [6] we briefly mention alter- 
native astrophysical explanations and also explain the 
prominent role of local sources. Finally, in Sec. Q we 
summarize the current status and outline some perspec- 
tives for the future. 



2. An excess... with respect to what? 

The claim of "anomaly" in the positron fraction has 
been initially presented by the PAMELA collabora- 
tion as well as in many interpretation papers as an ex- 
cess with respect to the no-reacceleration model pre- 
diction of [20], although with some caveat concern- 
ing the absence of error estimate. At face value, it is 
unclear to what extent one can consider the positron 
data as "anomalous", unless some firm statement on the 
positron fraction prediction can be made. In 12111 . the 
anomaly was reformulated as a spectral feature (namely 
a rising curve) rather than as an excess problem. 

Let us illustrate this point more quantitatively. The 
positron fraction can be written as 



f(E) = 



1 



1 



1+(<D,-/«V) 



l+*4 eV 



(1) 



where the fluxes O,- refer to the ones at the top of the 
atmosphere, and a simple fit to the data above 10 GeV 
(to minimize the charge-dependent solar modulation ef- 
fects) yields g = -0.38 + 0.06 for the original PAMELA 
analysis |3|], or g — -0.23 + 0.04 for the alternative 
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analysis presented in [4]. At the same time, a fit of the 
PAMELA electron data over the corresponding energy 
interval provides <& e - oc #-3.23±o.02 As a consequence, 
one infers an absolute positron spectrum which goes 
as OV oc £-2.85±0.06 ( reS pectively <D e + oc £-3.00±o.04^ 

Now, the proton flux measured by PAMELA itself in 
the similar energy -range is Op oc £- 2 - 82 M23I1 . and it is 
this hadronic CR spectrum which enters as source for 
the ISM positron one Q. In a pure diffusive propaga- 
tion framework, secondary /primary yields indicate that 
diffusion alone would steepen the spectrum by about 
6 - 0.4, see e.g. B2411 for more quantitative details. 
More importantly, TeV-leptons are subject to large en- 
ergy losses, steepening their spectral index significantly 
more: for example by (1 + 8)/2 « 0.7 in a two-zone 
diffusion model (where sources are embedded in a thin 
disk immersed in the thick halo f25lo or by ~ 1 in a 
leaky box model where radiative losses fully dominate 
the propagation. For illustration, the recent numerical 
analysis in 1 26] predicts a secondary positron spectrum 
roughly declining as <t> c + oc E~ 3 A 4- E~ 3 - 5 in the range 
of interest, clearly yielding a e + fraction inconsistent 
with the data. In view of some concerns about theoreti- 
cal assumptions entering these considerations raised for 
example in B2711 . it is worth noting that the only "the- 
oretical" argument is to assume that, as a result of the 
diffusion-loss process in the propagation through the 
ISM, there is no viable (i.e. consistent with known 
physics and cosmic ray astrophysics) mechanism which 
can steepen the spectral index of 10-100 GeV positrons 
with respect to their production by ~ 0.2 or less, i.e. by 
less than the purely diffusive steepening effect. Need- 
less to say, should such a kind of alternative solution be 
found, it would constitute a much greater breakthrough 
than the addition of primary sources @. Of course, the 
need for an additional component can be confirmed in 
more refined statistical analyses, see M M 51, but 
we believe that it is important to stress the simplicity 
and generality of the basic argument. 



'The recent indication by CREAM |2l and PAMELA (2J that 
Helium spectra are about ~ 0.1 harder than proton ones and that a 
hardening takes place in the 10 2 -r 10 3 GeV range at several hundreds 
of GeV does not change qualitatively the conclusion. 

2 Note that in the "nested leaky box model" the interpretation of 
the data is completely different, but one assumes as starting point that 
in fact there are fwo components of "secondary" CRs like Li, Be, B, 
e* and p (see 12811 for details). These are produced both in the ISM 
and inside high density "cocoons" surrounding sources, successively 
escaping in the ISM. We do not address here advantages and problems 
of this model, but note that in that case the need for a source contri- 
bution to "secondaries" is built-in as assumption. Of course, debating 
on the need of a source component is only meaningful in diffusion or 
leaky box models trying to fit the data with only ISM secondaries. 



Since we argued that the problem is inherent to the 
positron spectrum, present data strongly disfavour alter- 
native scenarios where the anomalous positron fraction 
is actually attributed to a non-conventional shape of the 
electron flux, e.g. due to the sum of a diffuse disk popu- 
lation and a local one as in [32]. While it appears natural 
that a homogeneous approximation for lepton propaga- 
tion breaks down at high -energies — even below the TeV 
scale — and likely future high-statistics data will be sen- 
sitive to inhomogeneity effects, quantitatively they ap- 
pear too small to explain the anomaly in present obser- 
vations (see also the recent analysis [33] suggesting that 
the inhomogeneities assumed in B32I1 are too extreme.) 
We come to the conclusion that the simplest and most 
conservative explanation of the observed phenomena is 
that a population of primary positron (or more likely of 
positron and electron) accelerators exist. 

This conclusion is gaining further support from re- 
cent Fermi measurements of the positron flux (i.e. sep- 
arating it from the total e~ + e + flux) via deflections 
in the Earth magnetic field, reported in [15], revealing 
a positron spectrum certainly harder than £ _3 -° in the 
range of interest. One component fits (i.e. only ISM 
production) of Fermi data fail to reproduce both the to- 
tal e~ + e + flux as well as the e + flux [34], providing 
additional support for: a) dismissing proton contamina- 
tion as main source of the effect; b) the need for an addi- 
tional electron-positron source to provide a satisfactory 
fit of the data. 

Once accepting that in order to explain the data one 
needs a primary source of positrons in addition to the 
secondary ones produced in the ISM, the pressing ques- 
tion becomes: what is (are) this (these) source(s)? Many 
answers have been proposed in the literature and we 
shall try to review them in the following, with an em- 
phasis on the likely astrophysical explanations. 

2.1. Another anomaly? 

Although way less discussed than its high-energy 
counterpart, by simple inspection of Fig.Q]one notices a 
discrepancy at low-energy between PAMELA data and 
previous experiments like AMS-01 and HEAT. In a cer- 
tain sense, this is even more puzzling: an experimen- 
tal agreement is found in an energy range which dis- 
agrees with the simplest expectations, while disagree- 
ment among experiments seems to appear where no new 
source effect is expected! 

More seriously, this trend might even raise the ques- 
tion if the entire measurement is plagued by some sys- 
tematic problem. It is not the purpose of the present 
review to discuss this issue in details. However, it suf- 
fices to say that it is generally accepted that no incon- 
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sistency is present at low energy, the effect being due 
to charge dependent solar modulation. Some investi- 
gations of this effect have been presented e.g. in ll35ll . 
and it has been extensively studied in l36ll which we 
address for further details. Since PAMELA has been 
taking data for half a solar cycle, it should be possible 
by now to see this trend to vary when time-binning is 
used, as well as an anti-correlation with a similar ef- 
fect in antiproton/proton data. It is clear that with more 
and more high-quality data available over long period of 
times, antimatter studies are becoming a more and more 
refined tool for heliospheric studies, see for example the 
forecasts for AMS-02 presented in l37ll . 

3. Is DM a good bet for explaining the CR lepton 
features? 

At first sight, an "excess" in the positron fraction 
is qualitatively expected in DM annihilation models, 
hence explaining the excitement generated by these data 
among particle physicists. Unfortunately, typical pre- 
dictions from WIMP dark matter annihilation fail to re- 
produce the data in three main aspects: 

i) First of all, the peculiar signature for an exotic ori- 
gin of the signal would rather be a spectral edge 
(more or less sharp, depending on the final state) 
after an initial rise; this drop back to the expected 
background at high energies is not observed, nei- 
ther in the e + fraction nor reflected as a dip in the 
e + + <?~ data[| 

ii) Insisting in attributing the rise to DM, the normal- 
ization of its contribution seems to be surprisingly 
large, compared with typical expectation for a S- 
wave annihilating thermal relic matching the ob- 
served dark matter abundance: i.e. in terms of an- 
nihilating cross sections, (crv) » (cv)th s-wave ~ 1 
pb. 

iii) Such a large yield of e ± pairs should imply large 
yields of high-energy particles in other channels; 
in particular one expects anomalies in antipro- 
tons, gamma-rays and possibly neutrinos, while no 
anomaly has been revealed yet in this respect. 

The first kind of problem can be overcome by pushing 
the mass of the DM nix to values sufficiently higher than 
the range explored by PAMELA, perhaps to the few Te V 



3 The so-called "ATIC peak" |6] has not been confirmed by later 
Fermi 

d,Q3] and HE SS data dUl, although no consensus exists in 
the community for an explanation of the disagreement. 



scale where the lepton fluxes appears to fall (see Fig. 
2). Since the differential spectrum from DM scales as 
m x 2 , this exacerbates the second problem and, as we 
shall comment upon, also affects the third one. In turn, 
the second problem can be cured in three qualitatively 
different ways: 

1. By attributing the high flux to a nearby clump of 
DM (as proposed for example in B38IO . which ap- 
pears however very unlikely in most structure for- 
mation scenarios and given existing gamma-ray 
constraints |39m . 

2. By accepting that (crv) » 1 pb and trying to de- 
rive that in a consistent model. The most popu- 
lar scenario is to invoke a Sommerfeld enhance- 
ment [40], i.e. the effect of non-perturbative en- 
hancement of the annihilation due to long-range — 
compared to the TeV~' scale — attractive forces: 
this would be active in the low-velocity v ^ 10" 3 c 
environment of our Galaxy, but would be inop- 
erative in the early universe, thus making the 
standard calculation of relic abundance still ap- 
proximately valid. In general, to accommodate 
such a mechanism one needs to introduce addi- 
tional light scales/particles with special coupling 
properties, at the expense of large fine-tunings. 
Alternatives include: i) non-thermal dark matter 
candidates — whose production does not happen 
by freeze-out, rather e.g. from decays of other 
species/condensates; hence (crv) can be large — 
ii) some Breit-Wigner resonant enhancement; iii) 
to abandon standard cosmological evolution alto- 
gether. Of course, these options still require some 
ad hoc tuning of parameters, since the relic abun- 
dance of the particle is decoupled from its indirect 
signatures. 

3. By considering a mechanism of production of the 
positrons which does not involve annihilations at 
all, rather DM decay. Unfortunately, the lifetime 
^ 10 26 s required to fit the data is not predicted 
in these scenarios, rather obtained a posteriori, al- 
though its scale might be justified in GUT models. 
The same "decoupling" from production affecting 
solutions of type 2 is of course plaguing these sce- 
narios. 

It is impossible to review here the different models 
proposed: a brief (and by now incomplete) review on 
the subject was attempted in 14 ill . In any case, the ab- 
sence of signals in channels other than the e ± and the 
type of spectrum observed pose serious constraints on 
the particle physics details of the models. Since the 
presentation of first diffuse data by Fermi, more and 
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more stringent limits have been derived from gamma 
rays, starting from the exploratory studies 11421 14311 to 
the more recent ones in [44, 45], which add to stringent 
bounds following from cosmic microwave background 
data (for a recent study see l46ll ). antiprotons IHtIi . or 
radio data 1481 14911 . Models that evade these constraints 
are quite contrived, see e.g. 51], and are unlikely 



to survive a closer scrutiny since they involve many 
approximations. A possibly incomplete list of these 
is: a) they do not provide a fit to diffuse gamma-ray 
data, i.e. they implicitly need astrophysical sources any- 
way, whose role would be to tighten the constraints; 
b) for heavy DM, they neglect effects related to elec- 
troweak bremmsthralung B52i l53l l54i 155 , 



, which 

have been shown to impact constraints (see for exam- 
ple 1 57, 58, 59, 6^]); c) when engineering "light medi- 
ators" to have a mass just below some hadron mass to 
forbid significant annihilation into gammas and antipro- 
tons, three-body final states are neglected, although in 
some concrete examples this has been proven to be a 
bad approximation 116 1[ 16211 : d) in the case where a light 
scalar <p with relatively large coupling gx to DM is the 
particle emitted in annihilation final states, a significant 
correction is expected to the two body channel. For ex- 
ample, the 3-<p emission with a soft <p should yield rel- 
atively large "infrared" log corrections to cross section, 
of the kind g x /{4n) loginix/m^). This in turn softens 
the spectrum and enhances the "boost" required to fit 
a given dataset of leptonic data at high-energy, besides 
making the model more vulnerable to low-energy CR 
constraints. 

Note that if one sticks to the conventional WIMP 
paradigm for DM, the present data by themselves are 
rather neutral to it, since typical expectations for DM 
signals in antimatter fall a couple of orders of magni- 
tude below the fluxes observed. We do not treat the 
DM explanations for CR leptons any further, but we 
note that the burst of activity following the data had 
at least two healthy effects on the community: on one 
hand, many reflections have been developed on hidden 
assumptions in the standard WIMP paradigm, and the 
spectrum of theoretical possibilities have been broad- 
ened. On the other hand, a reality check has been im- 
posed on the whole program of "indirect DM detec- 
tion"; this has revealed clearly that, facing the poorly 
known astrophysics, barring some exceptions no indi- 
rect evidence for DM can be robustly found outside 
a well-defined and overconstrained framework such as 
the vanilla WIMP paradigm. Current data are already 
limited by systematic uncertainties on the astrophysics, 
rather than by statistics (see e.g. the case of the diffuse 
emission ll45ll or Galactic Center J63ll one in gamma- 



rays.) It is fair to conclude that, until a better under- 
standing of the astrophysical sources is achieved, most 
antimatter signatures of WIMPs are far from robust and 
a "blind" search via multi-messenger data mining ap- 
pears poorly motivated. A new chapter might be even- 
tually opened in DM phenomenology if the search for 
CR traces could be formulated as an a priori problem, 
namely with input from a discovery at colliders and/or 
direct detection experiments underground. In that case, 
in principle correlated predictions among many chan- 
nels/directions/energies can be made which are unlikely 
to be mimicked by uncorrelated astrophysical sources. 



4. Pulsars (and Pulsar Wind Nebulae) 

The idea that pulsars might be associated with the 
production of cosmic ray electrons/positrons is quite 
old, see 16611 . and has been periodically recon- 

sidered in the past decade, see e.g. B671 16811 . The reason 
is that, being pulsars identified with fast rotating mag- 
netized neutron stars, a large electric field is induced 
which can extract electrons from the star surface: in fact 
these effects are so strong that a pulsar rotating with an- 
gular velocity Q is not living in vacuo, rather it is sur- 
rounded up to a distance known as light radius ri = c/Q 
by a comoving plasma configuration called "magneto- 
sphere" (see e.g. M Ez3, ffih •) 

In turn, the stripped electrons lose energy via curva- 
ture radiation (plus additional processes) while prop- 
agating far from the star along the magnetic field 
lines, and the emitted photons are so energetic that an 
electron-positron pair can be formed in the intense neu- 
tron star magnetic field. Numerous QED processes in- 
duce a multiplicative cascade which populates the mag- 
netosphere of the pulsar with pairs. The pairs produced 
in the magnetosphere, together with the Poynting flux 
emanating from the pulsar, form a relativistic magne- 
tized wind evolving in a rich environment: since the 
pulsar is born from the collapse of a massive star (core 
collapse supernova) it lies initially well within the ejecta 
of its progenitor, in turn surrounded by the supernova 
blast wave propagating in the ISM. When the cold, mag- 
netized relativistic wind launched by the star hits the 
non-relativistically expanding ejecta, a shock wave sys- 
tem forms in the impact: the outer one propagates in 
the ejecta, while a reverse shock propagates back to- 
wards the star. The latter is known as termination shock, 
where the wind is slowed down, its bulk energy dissi- 
pated and turned into that of a relativistic ally hot, mag- 
netized fluid, which then shines as a Pulsar Wind Nebula 
(PWN). 
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If we denote with B s the surface magnetic field and 
R s the radius of the neutron star, from Faraday's law 
one can naively estimate a potential drop of up to 
ip a>B s R 2 ~ 10 16 V available for accelerating e ± 
in the magnetosphere. Those pairs (which are more 
closely associated with the coherent pulsating radio and 
gamma emissions) suffer at least adiabatic energy losses 
while reaching the termination shock. Nonetheless, at 
the termination shock a relatively large fraction (a few 
tens of percent) of the wind bulk energy is converted 
into accelerated pairs, which then radiate a broad-band 
photon spectrum, extending from radio frequencies to 
multi-TeV gamma-rays, through synchrotron and In- 
verse Compton processes (see [72] for an overview.) 

It is unclear what the ultimate fate of these pairs is: 
while propagating far from the termination shock ad- 
vected together with the toroidal magnetic field lines, 
one expects them to be confined within the cavity cre- 
ated by the wind and to progressively lose energy. The 
fraction of them escaping in the ISM depends probably 
on astrophysical details. 

Given this broad picture, can one attribute the sig- 
nal inferred from the charged lepton data to the PWN 
contribution? There are several elements which suggest 
that a positive answer is possible if not likely, which we 
briefly address below. 

4.1. Energy budget 

A clear upper limit to the total energy into e ± is given 
by the pulsar spin-down power integrated over time, i.e. 
the energy loss corresponding to the slowing rate of ro- 
tation which in the standard theory provides virtually 
all of pulsar energy. In sufficiently general form, the de- 
crease of the rotation frequency £!(= 2n/P,P being the 
period) can be written 



Q 



-a£l\ 



(2) 



yielding a spin-down luminosity (/ being the momen- 
tum of inertia) 



L= \E\ = /£l|£l| = a/£l" +1 
The solution of Eq. (0 yields 
£1 



(3) 



£1(0 = 



£1" +I 

£{t) = al (4) 



where we introduced the characteristic timescale 



to = [a(n - mVr 1 



(5) 



The special case of magnetic dipole corresponds to 
n = 3 and a = 5 B 2 S R A J{% M s c 3 ), where M s the mass of 



the neutron star. By integrating the luminosity over time 
one gets 



-ml 

2 



2.2 x 10 4 



(6) 



I AM, 



W R s \ 2 /£lp\ 
jllOkm/ \Hzj 



erg, 



which amounts typically to E, , 10 49 erg or more 
mostly converted into a magnetized, relativistic wind in 
a time comparable with tq. Since one expects a rate 
Rcc °f about 2 core-collapse SN per century in our 
Galaxy (following for example from 26 Al gamma-ray 
data B73I1 ). the maximum luminosity injected is of the 
order of 

Rcc E tol 



m,* ~ 6.3 x 10 3y ergs~ 



(7) 



2 century^ 1 10 49 erg 

to be compared to what needed to fit the data which is 
about one to two orders of magnitude lower, see e.g. |74 , 



4.2. Spectral shape 

Observationally, the spectrum of radiation from 
several PWNe requires a spectrum of electrons and 
positrons which has a broken power-law shape: the 
break happens at Lorentz factor of 10 5 4- 10 6 — hence 
at the (9(100) GeV scale — with a very hard power-law 
below the break, of the kind ^ E~ l 5 , while significantly 
softer (softer than E~ 2 ) above the break (see e.g. 11721 
p76ll ). It is important to emphasize that this is only 
based on data (mostly radio-to-X ray data) and does not 
rely on theoretical arguments. The hard spectral index 
matches what needed to fit the data 174l l75l l77l 17811 . 

So, on general grounds, one can conclude that: i) pul- 
sars have more than enough rotational energy to explain 
the normalization of the cosmic ray lepton features, 
ii) Evidently, a significant part of this energy is con- 
verted into kinetic energy of lepton pairs at the termi- 
nation shock of the PWNe, with the right spectral prop- 
erties to match the high energy positron fraction as well 
as e~ + e + energy spectra. Also, it is worth adding that 
observations from Fermi-LAT instrument show clearly 
that pulsars (with or without the identification of PWNe) 
are the most abundant population of Galactic objects 
to shine in the GeV gamma-ray band [79]. Including 
their contribution to any account of the leptonic CR 
sources is thus a requirement for any realistic descrip- 
tion of Galactic CR sources. Reversing the argument, 
it is somewhat surprising that until recently — when suf- 
ficiently precise leptonic data in the 100 - 1000 GeV 
have become available — one could fit the data without 
necessarily including them! 
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4.3. Theoretical open problems 

Of course, the above arguments do not settle the 
many theoretical issues that still remain in the interpre- 
tation of these objects, for details see for example [80]. 
What PWNe observations seem to imply is a relativistic 
wind mostly made of electron-positron pairs (so mostly 
carrying particle kinetic energy), although it originates 
as a Poynting flux-dominated outflow in the magneto- 
sphere close to the pulsar. How does this conversion 
take place? In particular, how and where does the large 
multiplicity of pairs implied by observations originate? 
The broken power-law spectrum is also very peculiar: 
even for the high-energy part which is steeper than E" 1 , 
in apparent agreement with expectations from diffusive 
shock acceleration at a relativistic shock, the highly rel- 
ativistic quasi-perpendicular nature of the shock should 
make acceleration inefficient, see e.g. 181]. The very 
hard low-energy spectrum is even more puzzling, since 
it implies a situation where most non-thermal particles 
lie at low-energy, but the energy is carried to a large ex- 
tent by the few particles close to the high-energy break. 
To name but one of the associated puzzles: why the (ex- 
pected) bulk of thermal low-energy particles acting as 
reservoir for the accelerated ones has never been identi- 
fied in PWN? 

Still, several interesting ideas have been put forward 
which address at least partially some of the issues. A 
cyclotron absorption acceleration mechanism ifs^h has 
been proven in simulations ll83ll to be able to produce 
the desired spectra and efficiencies if a substantial frac- 
tion of the pulsar wind energy is carried by protons or 
other ions, namely if they are energetically dominant 
even if their number is negligible with respect to lep- 
tons 0. Interestingly, it could also present preferential 
acceleration of positrons thanks to helicity-matching of 
the waves onto which they scatter, which are generated 
by the (positively charged) protons. But the required 
high Lorentz factor of the wind, although to some ex- 
tent theoretically motivated, seems at odd with recent 
simulations as e.g. in 117611 . 

Another conceivable mechanism invokes magnetic 
reconnection of the striped magnetic field of the wind at 
low latitude around the pulsar rotational equator, which 
at the termination shock would convert into particle ki- 
netic energy 18411 . The shape of the particle spectrum is 
however dependent on poorly understood details and it 
is difficult to confirm it observationally at present. Yet, 
recent two- and three-dimensional particle-in-cell sim- 
ulations seem to support this kind of scenario [85]. 



4 This is of course conceptually possible, given the large rest mass 
ratio m„ jm e 1836. 



It is worth emphasizing once again that, although it 
is mandatory to solve these issues in order to put the 
theory of PWN acceleration and non-thermal emission 
on more solid ground, to some extent these problems 
are less directly linked to the contribution that PWNe 
can give to CR leptons seen at Earth. Since the spec- 
trum of PWN is relatively well-known, more crucial 
to answering cosmic -ray questions are issues rarely ad- 
dressed such as the escape probability of the accelerated 
pairs into the ISM, simply because they are way less 
relevant for "photon astronomy". In that respect, an in- 
teresting proposal has been done in Ref. 18611 . where it 
is argued that the so-called Pulsar Bow Shock Nebulae 
(i.e. high-velocity pulsars which have escaped their host 
SNR) might release their high-energy leptons directly in 
the ISM, a few tens of thousand years after their birth. In 
general, since the efficiency of conversion of spin-down 
luminosity into CR leptons is quite high, the effective ef- 
ficiency of 0{ 1 %) would follow mostly from the smaller 
residual energy fraction available after the escape from 
the SNR, rather than from a relatively inefficient conver- 
sion. This argument is also at ease with the non-dipolar 
spin-down laws — 2 < n < 3 in Eq. (f2| — found in the 
few pulsars for which the second derivative Q has been 
measured (n - -Q.Q./Q. 2 ). 

5. Production in SNRs 

Usually, one neglects the production of "secondaries" 
within SNRs with the argument that the probability of 
interaction for a primary species like protons is small 
over the lifetime of the accelerator, and in any case 
smaller than the probability of interaction during the 
diffusion time in the ISM. For example, for a typical 
inelastic cross section of 30 mb one finds an interac- 
tion probability of about 10~ 2 5 in a typical ISM density 
of 1 particle crrT 3 and for a representative lifetime of 
a source of 10 5 yrs. This has to be compared with a 
probability more than one order of magnitude larger to 
undergo the same process during the propagation in the 
diffusive halo, which takes up to 0(1 7 ) yr, although in 
a medium whose density is in average much smaller. 

However, this argument does not take into account 
that the spectrum resulting from the ISM production 
is steeper than the source one due to diffusion effects 
(roughly E~ 2J vs E~ 22 ), because it depends on the 
steady state spectrum of primaries rather than their in- 
jection one. Hence, one expects that at sufficiently high 
energy the source production terms can make a signif- 
icant contribution to the overall secondary spectrum, 
as already noted in [87] for gamma-rays and more re- 
cently discussed for antiprotons ll88ll . Naively, an ap- 



7 



proximately flat secondary/primary ratio is expected as 
a result. 

Recently, it was noted in [89] that an even more pecu- 
liar behaviour may follow when realizing that the sec- 
ondary production takes place in the same region where 
cosmic rays are being accelerated. With the reason- 
able assumption that the diffusion coefficient D of rela- 
tivistic particles near the accelerating SNR shock grows 
with energy, the region from which particles can re- 
turn to the shock and undergo further acceleration gets 
larger and larger with energy. If, at the same time, the 
maximum energy attainable E m;lx is not limited sim- 
ply by the naive diffusion timescale t&e ~ D/v 2 (v 
being the shock velocity scale), one may obtain grow- 
ing secondary/primary ratios, which can fit the current 
data. Note that a non-trivial behaviour of E nmx is needed 
anyway to allow particles to reach knee energies via 
diffusive shock acceleration, since the naive estimate 
from typical ISM values for the diffusion coefficient 
would yield much smaller energies, as noted a long time 
ago B90I1 ! A lot of attention has been paid lately on the 
proposal that this problem could be overcome via a non- 
resonant amplification of magnetic field by CRs them- 
selves [91], which seems to be confirmed (although not 
unambiguously) by the relatively thin X-ray rims seen 
in several SNRs, e.g. Ir92ll . In general, a firm predic- 
tion is hampered by the very uncertain evolution of the 
magnetic field (amplification, damping, etc.) especially 
at the late stage of SNR evolution which are thought to 
contribute the most to "low-energy" (i.e. sub-TeV) CRs. 
Additionally, non-trivial time and space evolution ef- 
fects are expected (both for the plasma dynamics around 
the shock front as well as the background "target" mat- 
ter), as for example recently tried to model phenomeno- 
logically in 19311 . Finally, the details of the injection 
and escape of CRs in/out of SNRs are yet unclear. So, 
the steady state effective description attempted in rfsoh 
is clearly a simplification of the picture only meant as 
effective model. Not surprisingly, time-dependent sim- 
ulations where most of these simplifications are retained 
but E max is fixed to its naive expectation find "only" 
a flat contribution to the Secondary /Primary ratio [94], 
while they confirm the rising behaviour found in [89] if 
£max is left free to assume larger values. Fortunately, 
this scenario appears to be testable independently of the 
fact that a firm theory is not yet available. In fact, if £ mal 
(or an equivalent quantity) is treated as a free parame- 
ter to fit e.g. the positron data, a clear prediction is ob- 
tained for associated observables like the antiproton flux 
or antiproton/proton ratio [88] or secondary/primary nu- 
clei 119511 . While the former predicts a rise at energies 
just beyond those currently explored, the nuclear ratios 



show a more ambigous picture, with some species like 
the titanium-to-iron ratio from ATIC showing evidence 
for a rise, and others like the B/C data from other exper- 
iments favouring a decreasing function. Most likely this 
issue will be settled in the near future, especially thanks 
to AMS-02 data (see for example the forecasts in [96]). 



Of course, the "textbook" model of SNR sources of 
CR via a strong, non-relativistic, planar shock is a "ze- 
roth order" approximation of what happens in nature. 
Not all SNe are identical, and significant qualitative dif- 
ferences may exist between relatively light and rela- 
tively heavy progenitors. Massive stars like Red Su- 
pergiants (15M < M < 25M Q ) and Blue Supergiant 
(or "Wolf-Rayet") stars (M > 25M Q ) explode into their 
stellar wind, which is magnetized, relatively dense and 
helium and metal-enriched from exposing the deeper 
layers of the star through mass ejections. The magnetic 
field topology is expected to be radial in a polar cap, 
and tangential over most of the remaining solid angle. 
From the different conjectured forms of the diffusion 
tensor and acceleration mechanism in these two envi- 
ronments, it was argued that a harder (E~ 2 ) component 
from the polar cap would eventually emerge over the 
steeper (£T 7 / 3 ) component which is numerically domi- 
nant at lower energies 19711 . By adjusting the normal- 
ization parameters, it was shown in 19811 that this al- 
lows one to fit leptonic data, and argued in lr99Tl that 
this is consistent with spectral breaks recently found in 
hadronic CR data II221 12311 . Qualitatively, distinguish- 
ing between this proposal (two alternative acceleration 
mechanisms within the same source, which in turn is 
the high-mass end of the SN sequence) and the above 
conjecture ("universal" albeit subtle features of DSA in 
SNR) might be difficult. One possible feature is that 
in the Wolf-Rayet scenario there is a potential contribu- 
tion to positrons from f} + unstable isotopes produced via 
photo-dissociation and spallation of the (abundant) nu- 
cle i (fur ther considerations on this subject can be found 
in llOOIO . Hence the normalization of rises in positrons 
and antiprotons or secondary/primary nuclei might be 
less sharply linked one another. Of course, a better un- 
derstanding of additional observables (isotopic anoma- 
lies, primary spectral shapes, etc.) may also help clari- 
fying the detailed scenario. Note that additional sources 
of e ± discussed in the literature involve the interac- 
tion with dense background photon de nsities possibly 
present near the accelerator llOll 1 1021 with more or 
less reasonable parameters, although it is unclear how 
the resulting leptons might escape unaffected from the 
sources while overcoming the strong radiative cooling. 
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6. Alternative explanations and contributions from 
local objects 



The well known Shakespeare's quote There are more 
things in heaven and earth, Horatio, Than are dreamt 
of in your philosophy [Hamlet Act 1, scene 5] cannot 
be more appropriate to describe the zoo of objects char- 
acterizing non-thermal astrophysics. It comes with no 
surprise that most of them have been associated to lep- 
tonic CR spectra. 

For example, some accretion-powered X-ray binaries 
made of a black hole accreting from a stellar compan- 
ion possess a relativistic jet and act as "micro-quasars", 
small galactic counterparts of their huge cosmological 
analogues. These objects are known to acc elera t e lep - 
tons at least to GeV energies during flares [103, 104]. 
According to some optimistic estimates, they may con- 
tribute to a sub-leadi ng ye t non -negligible level to the 
Galactic CR budget Jiol [ToS. In Ref. JTo3] such 
an object (or alternatively a Gamma-ray burst (GRB) 
with the right timing/distance) is mentioned as possible 
source of the anomalies. 

Similarly, in jl08|] it was speculated that white dwarf 
(WD) "pulsars'Q could potentially dominate the TeV e ± 
window. Although a number of such objects should ex- 
ist, their "detection" as non-thermal sources is still spec- 
ulative, not to speak of their frequencies and efficie ncies 
as accelerators. Yet another option investigated in II 1 0911 
is that positrons originate from collisions of hadronic 
CRs from a SNR in a nearby dense cloud environment 
(see also 1 1 id 111 1 for similar consid eratio ns devel- 
oped more than two decades ago). In ill 12ll . the role 
of magnetars (ultramagnetized neutron stars) is consid- 
ered, in addition to "ordinary" pulsars. Needless to say, 
the list is incomplete! 

From a more general perspective, why such a pro- 
liferation of source candidates? There are two rea- 
sons for that. First of all, it has been noted since 



in the relevant £ -range is due to radiative losses, and 
can be estimated in the Thomson limit as 



some time 11 131, ll 1411 that at higher and higher ener- 
gies the contribution of local leptonic CR sources be- 
comes more and more important, making the approxi- 
mation of continuous source distribution inappropriate 
(in particular above (9(100) GeV or so). Turning the ar- 
gument around, the TeV range for leptons has a large 
potential for astrophysical diagnostics, since it may re- 
vea l the contribution of single sources, as discussed e.g. 
in ill 15L This can be quickly understood in physical 
terms: the cooling time r coo i = E / (-dE / dt) for leptons 
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where u tot is the total energy density in low-energy pho- 
tons (CMB, infrared, optical . . . ) plus the energy den- 
sity in fi-field, ub = B 2 /8n. This implies an effective 
"diffusion horizon" 7?diff of the order of 



(9) 



i.e. sub-kpc distances for TeV electron sources. Sec- 
ond, the higher in energy the single source starts to con- 
tribute, the less demanding the energetic requirement 
are, although the requirement on E m . dx and the accelera- 
tion mechanism tighten. 

Given the above arguments, provided one finds a 
mechanism to accelerate leptons up to sufficiently high 
energies and with a sufficiently hard spectrum (and 
make them escape!), one can always envisage arrang- 
ing the distance, time, or normalization of the injec- 
tion in order to fit the lepton CR data. Of course, 
when considering known objects whose distance, age 
and rough energetics is known, as is the case for pul- 
sars, the exercise results significantly more constrained. 
Many phenomenological attempts have shown that is 
possible to fit the lepton data up to the highest ener- 
gies, by addingthe contributions of a discrete distribu- 
tion of pulsars Q which can be taken fro m (admittedly 
incomplete) pulsar catalogues ill 16l ll 17[ 1 1 1 Stl like the 
ATNF |119], some theoretical models for the distribu- 
tion of sources in the Galaxy (see e.g. 1 120l 12lll ). or a 
few prominent nearby so urces, like Geminga or Mono- 
gem in fzllizl EH DTI. 

Not surprisingly, when fitting 
the data with one or few pulsars a significantly larger ef- 
ficiency (typically a few tens of percent) of luminosity 
conversion into pairs is needed, which is probably also 
indicative of catalogue incompleteness — at very least 
for geometric reasons in the case of the pulsar, which 
can only be discovered when showing the right "beam- 
ing" towards the Earth. 

On the contrary, when invoking putative objects, 
like a past GRB, the plausibility of the contribution 
should rather be formulated in terms of probability 
that the "right" distance and timing conditions a re re- 
alized. For example the proposal formulated in 1107 1 



I.e., a rapidly spinning WD thought to be formed by a merger of 
two ordinary WDs, since the observed WDs are usually slow rotators. 



In some cases a lower cut on age is introduced for discarding 
sources where (presumably) leptons are still confined in the source. 
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(see also B123IP requires a GRB happening in our Galac- 
tic neighborhood in the last ~ 10 s -r 10 6 yr, whose 
chance probability is probably no higher than 0(1%) 
(even forgetting about the evidence that the Milky Way 
may be too metal-rich to have hosted bright GRBs re- 
cently II 2411 .) An additional aspect is that the cosmic 
ray spectra produced in (and escaping from) GRBs, mi- 
croquasars or more exotic sources are way less con- 
strained, which makes the case for invoking these ob- 
jects less compelling from the phenomenological point 
of view, while of course giving more freedom from the 
astrophysical model-building perspective. 

All in all, while it is not challenging to fit the positron 
fraction data with a continuous source distribution, see 
e.g. 1174 17511 . it is very likely that the high-energy part 
of the lepton CR flux reflects to large extent a few 
prominent sources; recently, this has been nicely illus- 
trated in II 1 2511 . with the help of semi-an alytic al meth- 
ods and montecarlo simulations (see also jl26l0 . A sort 
of "galactic variance" allows one to consider the range 
E >TeV for cosmic ray leptons as a window to "sin- 
gle source contributions": they can certainly be of the 
same type of those likely to contribute at low energy 
(like SNRs and PWNe) but possibly — although with 
smaller probabilities — of a special type. Clarifying this 
kind of astrophysical questions is the main purpose of a 
dedicated experiment, the "Calorimetric Electron Tele- 
scope" (CALET) 112711 which is designed to observe 
electrons up to 20 TeV (and nuclei up to 1000 TeV) 
and is expected to be placed at the Japanese Experiment 
Module at the International Space Station in the current 
decade. 



7. Perspectives for the future and conclusions 

Compared with the relatively slow progress of the 
previous decades, the last few years have seen a fast 
pace of new cosmic ray data. The qualitative improve- 
ment has been particularly spectacular for leptonic CRs. 
Theoretically, it is challenging if not impossible to ex- 
plain the data without assuming a primary source of 
positrons (or e ± pairs). Despite the early enthusiasm 
about possible dark matter interpretations, it proves very 
difficult to reconcile all the cosmic ray data with such an 
exotic explanation, which (besides having to face the- 
oretical challenges) is at present even phenomenologi- 
cally disfavoured. Perhaps more importantly, since at 
least a couple of decades it has been realized that as- 
trophysical objects exist which can produce the kind of 
features presently observed with high statistics, and al- 
ready hinted to by some past data: even limiting one- 



self to astrophysical candidates, numerous options have 
been suggested. 

Here we tried to provide an overview of the situa- 
tion following a basic criterion: When entering an "un- 
charted territory" of nature, just like at the beginning 
of the scientific era, in order to progress in our under- 
standing it is crucial to make use of Ockham's razor to 
limit the number of hypotheses, rather than adopting the 
strategy to explore all logical possibilities, however tiny 
their probabilities to be realized in nature are. In fact, it 
turns out to be already challenging to explore the con- 
sequences of playing with the tools at hand, namely to 
study more in depth the predicted signals deriving from 
objects known to exist (and actually being quite com- 
mon) in the Galaxy. 

This exercise singles out two leading classes of 
sources, SNRs and PWNe, whose energy budget is 
large enough to account for the data and are seen to 
shine in the non-thermal sky from radio-waves to GeV 
and TeV bands. The first problem with SNRs — which 
are nonetheless the prime candidate as hadronic CR 
sources — is that positrons are not present in the nonrela- 
tivistic "background plasma" into which the shock front 
propagates, so they must be secondaries. Inelastic colli- 
sions of protons/ions with the background medium may 
be enough to provide the sufficient number of positrons, 
and in some scenarios positrons may also be byproduct 
of unstable nuclei or reactions onto photons. A more 
challenging question is how to get a rising positron frac- 
tion rather than a flat ratio at most. In general, some ex- 
tra ingredient in the acceleration mechanism is needed, 
on which different models differ; also, all models pre- 
dict correlated "rising" signatures in other byproducts 
like antiproton/proton spectrum or secondary/primary 
nuclei. The good news is that forthcoming data from 
AMS-02 should b e abl e to reveal these signatures and 
test the basic idea jl28ll . 

Interestingly, even if SNRs turn out not to be the 
main contributor to lepton anomalies, it appears rea- 
sonable to expect some (in principle measurable) con- 
tribution to leptonic and hadronic secondaries at high 
energy. An important lesson is thus that traditional in- 
dicators of DM (like excess of antiprotons) or propa- 
gation diagnostics (as both antiprotons and B/C) may 
have a less trivial energy behaviour than often thought. 
In a sense, this was already knowrQ, but somehow 
forgotten in the last one or two decades when, per- 
haps in connection with the development of numerical 



7 For example, similar considerations had beed done some time 
ago in t he c ontext of cosmic ray acceleration in the interstellar 
medium ll29ll. 
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tools for CR propagation (such as GALPROP 113; 
USINE HH, DRAGON iBl ). more attention has 
been devoted to the uncertainties introduced by the 
"propagation" parameters than those entering the accel- 
erating sources. Fortunately, more and more attention 
is paid to these issues, as well as to the impact that 
the primary spectral shape at high en ergy has on sec- 
ondaries HI EH EH [l35L fl36iTl37ll . 

On the other hand, pulsars (or more exactly PWNe), 
although traditionally neglected in the CR budget, ap- 
pear to have both the right energy and the hard spectral 
shape for leptons at the termination shock to naturally 
explain the lepton CR anomalies. Yet, the theoretical 
understanding of acceleration in these objects is less ad- 
vanced, and some details (as the escape probability in 
the ISM) appear important to assess quantitatively their 
contribution to interstellar CR spectra. If they provide 
the main contribution to the lepton anomalies, it is un- 
likely that future charged CR data will answer the many 
obscure points. More probably, progress will come 
from improvements into the modeling of their global 



emission spectrum from radio to gamma band B761I138 I 



as well as from simulations of the microphysics playing 
a role in these objects, as e.g. in 118 lll84[ 18511 . 

Whatever the sources of the lepton anomalies are, a 
dominant contribution of one or few local sources be- 
comes more and more likely at higher and higher en- 
ergies. This is a generic prediction of virtually any as- 
trophysical model which can in principle be confirmed 
in two ways: i) via high-statistics measurements of the 
spectrum, as those expected from CALET B127I1 . which 
should show "fine structure" spectral features due to 
inhomogeneity and stochasticity of sources, ii) Via a 
relatively large anisotropy (roughly at the 0.1% level 
at (9(100) GeV), which is expected from n earby ob- 
jects [139]; current upper limits from Fermi jl40| ap- 
pear close to the most optimistic predictions, see II 14 lL 
1 14211 . Some perspec tives for AMS-02 in this channel 
have been studied in II 14311 . 

Let us end this review with a word of caution and a 
subjective point of view on forthcoming progress. As 
usual, one cannot but stress on the importance of ad- 
ditional and high quality data. It is even possible that 
current or future experiments might completely revo- 
lutionize our understanding of the field, as for exam- 
ple could follow from the discovery of anti-helium nu- 
clei in AMS-02 or antideuterons (perhaps from DM) in 
GAPS 1 14411 . Even lacking such major breakthroughs, 
it is guaranteed that data from AMS-02 alone will 
narrow down uncertainties in the parameters (like pri- 
mary spectra at high energy, diffusion coefficient, halo 
heigth, etc.) entering theoretical predictions of the sec- 



ondaries B128I1 . However, it is already clear with the 
current data that going beyond the zeroth-order descrip- 
tion of CR datasets by introducing multiple sources un- 
avoidably brings huge degeneracies due to the large 
numb er of unconstrained parameters, see for example 
1 12111 . To some extent, going beyond the overall dis- 
tinction between classes of sources (like PWN vs SNR) 
to the level of "inverting" the problem (i.e. finding ac- 
tual sources given the spectrum) might make little sense, 
given the unavoidable incompleteness of high energy 
astrophysical catalogues and the fact that, due to diffu- 
sive propagation, some sources contributing now to CR 
data might not be visible anymore in other bands: the 
high-energy universe is in fact time -dependent. Thus, 
it is worth remembering that knowledge in this field 
does not automatically follows from the collection of 
more data: an example is provided by the puzzling ob- 
servations of (hadro nic) CR multi-TeV anisotropy from 
many experiments fill EH EH EH EH Em for 
which no convincing explanation has been found, yet. 
To put that in the perspective of our original problem 
of lepton diagnostics, detecting anisotropy might not be 
enough to make source identification, since other effects 
(like inhomogeneous diffusion coefficient, anisotropic 
diffusion, local bubbles, etc.) may even mask or alter 
the source effect. Perhaps energy-dependent anisotropy 
studies might be needed for diagnostics, but these ap- 
pear even more-challenging. If compared with the ef- 
forts gone into modelling cosmic ray propagation or as- 
trophysical models for photon astronomy, the theoreti- 
cal understanding of high energy sources of the cosmic 
rays eventually detected at the Earth (i.e. not remotely!) 
is much more rudimentary, also because of the highly 
indirect nature of the problem. It is auspicable that 
the excitement brought in the astroparticle community 
by the recent data will also stimulate a new generation 
of scientists to devote themselves to theoretical aspects 
of acceleration and escape of CRs from astrophysical 
sources, since firm progress in solving these fascinat- 
ing questions will certainly require additional ideas and 
tools and a critical re-examination of the many simpli- 
fying assumptions underlying present scenarios. 
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